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SUMMARY 


During the past few years, the direct-inverse technique has been de- 
veloped into a numerical method, called TRANDEL, that is suitable for the 
analysis and design of subsonic and transonic airfoils and for the evaluation 
of design concepts. This paper provides a general description of the method, 
demonstrates its application to a design-analysis type of problem, and 
finally, discusses a new usage of the method for the low speed high lift case. 


INTRODUCTION 


The basic concept of the present method (refs. 1-4) is to have a tech- 
nique which can be used In either the direct (analysis) mode In which the 
airfoil shape is prescribed and the flowfleld and surface pressures are de- 
termined, or in the direct-inverse (design) mode in which an initial nose 
shape is given along with the pressure distribution on the remainder of the 
airfoil. In the latter case, the flowfleld and actual airfoil shape are 
computed . 

The resultant computer program, called TRANDES, (ref. 5), has several 
unique features. In order to achieve accuracy, the method utilizes the full 
Invlscld potential flow equation; and in order to remain simple, it solves 
the problem in a stretched Cartesian grid system that maps the infinite 
physical plane to a rectangular computational box. Further, to avoid at 
supersonic points difficulties associated with nonalignment of the coordinates 
and the flowfleld, a rotated finite difference scheme is used in the solution; 
and the resulting transformed finite difference equations are solved itera- 
tively by column relaxation sweeping from upstream to downstream. Finally, 
the method does Include the effects of weak viscous Interaction. The basic 
idea in the design case is to treat the airfoil determined by the Inverse 
method as the displacement surface and to subtract from it a displacement 
thickness determined by a Nash-Macdonald (ref. 6) turbulent boundary layer 
computation in order to obtain the actual airfoil coordinates. It should be 
noted that the. present program determines the airfoil shape simultaneously 


*Supported by NASA Grant NSG 1174 

**Presently with the Boeing Co., Seattle, WA 


55 


with the flowfield relaxation solution. For the analysis case, the approach 
is to calculate a boundary layer displacement thickness and to use it to 
correct the location of the displacement surface (i.e. airfoil ordinate plus 
6*) • The flowfield is then solved iteratively with the displacement surface 
being updated every ten relaxation cycles. 

Now, in the design mode, the shape of the nose region (typically 6-10% 
chord) is specified and a pressure distribution is prescribed over the re- 
mainder of the airfoil. Thus, the appropriate airfoil boundary condition in 
the direct region near the leading edge is the surface tangency requirement 
and in the inverse region, where the pressure is specified, it is essentially 
the specification of the derivative of the perturbation potential in the x- 
direction. In order to satisfy these at the airfoil boundary, which in 
general will not coincide with the Cartesian grid points, the derivatives in 
the boundary conditions are expanded as two term Taylor series about a dummy 
point inside the airfoil. The derivatives in these series are then written 
in finite difference form using second order formulas for all first 
derivatives and at least first order ones for higaer derivatives. In the 
direct region, central differences are used for :c-derivatives and forward 
(on the upper surface) for the y-derivatives. However, to prevent numerical 
Instability, the inverse region uses a second order backward difference 
formula for the first term of the Taylor series representing the x-derivative. 
For details concerning the finite difference formulation, boundary conditions, 
etc., see references 2-5, and 7. 

Currently, the program can be easily used to design or analyze an air- 
foil at a specific flight condition. Figure 1 shows such a result for the 
design of a low lift airfoil having a sonic rooftop followed by a linear 
recovery. Supercritical shockless airfoils have also been designed and 
examples are presented in reference 7, 

Now any numerical method needs to have its accuracy verified. As a 
result, comparisons with other methods (refs. 8,9) have been conducted for 
airfoils ranging from 4 to 16% thick and Mach numbers up to 0.85. Some of 
these results are presented in references 2, 4, and 10; and, in general, 
the agreement is excellent. In addition, comparisons have been made with 
transonic experimental data obtained at the NAE (Ottawa) (ref. 7) and at the 
Ohio State University. A typical comparison with the OSU data is shown foj 
the GA(W)-2 airfoil in figure 2. In general, the agreement for C_, C , 
and is quite acceptable. ^ ^ 

Symbol definitions are given in an appendix. 

APPLICATION TO DESIGN 


To demonstiate the usefulness of the present combined design-analysis 
method, consider rhe problem of designing an airfoil having a rooftop plateau 
followed by a Stratford recovery (ref. 11). (It should be noted that a Strat- 
ford recovery is used here only for example purposes. While such a pressure 
distribution has the advantage that it maintains the boundary layer at a 
constant margin from separation, it also hi_s several disadvantages which will 


the addition, two design concepts will be compared. In 

e first case, the airfoil will be designed at a low with a sonic unoer 
rooftop pressure plateau; while in the second ^ case, the LrfSil 

win ^h 3 supercritical plateau at a moderate C, . The airfoils 

will then be a^lyzed and compared at conditions other than ^ the desien 
point to see which one has the better characteristics. 

0 35 rooftop design, which had a target C of 

0.35, are sho|m in figure 3. The solid line is the design C dlsLl- ^ 

iTtZ' I accounting for the^ boundLy layer 

r!c T reverse curvature on the upper surface Ld the 

thinning. This behavior is typical of airfoils employing 

?ura?frf frequently if: leads to shapes that are struc- 

turally too thin in the vicinity of the trailing edge. 

of ^h!°T obtained? Since the method requires the specification 

a region, the nose shape can be used to control the trailing 

edge thic^ess. This approach is demonstrated in figure A and works easiest 
if an analytical nose shape controlled by a single parameter is utilized. 

Here the nose shapes are those associated with NACA OOXX airfoils. Notice 

increases, the trailing edge thickens. In ad- 
should be noted that all the results in figure A were obtained 
using the same Cp boundary conditions and that each case is completely 
Independent. Thus, obvious physical unrealities such as trailing edge 

my*be"^used final design shape. Also, asymmetric nose shapes 

specifically in figure 5, there is a 
between ttis values of the nose shape parameters and the 
illng edge ordinates . Usually the variation is essentially linear; and 
us, a 2 r o taiiiing the results for two cases, the desired trailing edge 
ordinates can be obtained on the third try. * * 

Sometimes in an aft-camber design case if the computed upper and lower 
surfaces are not almost parallel near the trailing edge, the flow in that 
vicinity will deviate from the desired pressure distribution and try to 

and^Jo^rh 'tt separation. Usually by slightly adjusting back 

nd forth the starring point of the upper surface recovery, a location can be 
oun w c will yield acceptable trailing edge slopes and pressures. This 

IZll TrT. iterations and some adjusLent in nose 
shape, but it is normally not difficult. 

Figure 6 shows ti.“ design pressure distribution and resultant airfoil 
shape obtained when the .irfoll was designed for a higher C and with a Mach 
1.1 supercritical rjoftop. Since this airfoil was designed^ at higher lift 

tJV “ surface pressure distribution had a larger aft bucket. 

This case demonstrates the disadvantage of using a Stratford recovery on a 

of*80% chord^"*’^*^^'^ airfoil in that the airfoil thickness is less than 2% aft 
Now an essential requirement of a design method, from an engineering 
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standpoint, Is that when the designed airfoil Is analyzed at the design 
condition, the computed Cp distribution should agree with the Cp distribution 
used for the design. Such analysis results are shown In figures 1, 3, and 6 
and were obtained using the present, method with viscous Interaction Included. 
As can be seen, the agreement Is excellent. It is believed that these com- 
parisons verify the engineering consistency of the method and, since the 
analysis results usually agree with experimental data, that the airfoils 
designed by this method should perform as predicted. 


APPLICATION TO ANALYSIS 


One of the difficulties associated with using a Cartesian grid for an 
analysis computation Is that such a grid does not place a large number ot 
points near the leading edge. Thus, the wave drag coefficient, which Is 
determined by integration of the pressure coefficient, has an Inherent error. 
Previous studies have determined this error Is consistent and primarily a 
function of airfoil shape and grid spacing; and thus, correction factors 
can be determined from calculations at subcritical speeds where the wave 
drag should be zero. Unfortunately, the correction procedure suggested in 
reference 5 may be partially In error; and while research is continuing, 
the results presented In this section have bean obtained using the following 
technique. 


At a Mach number, which the flow is entirely subcritical, 

determine for each angle of attack, a, the normal, and axial, 

Ca , , force coefficients. Then find the drag correction from 

Asub,a 

“dw ^ -'a , ^ 

sub, a sub,u sub, a 


Next find the correction at the ueslred supercritical Mach .number, M , by 


L-nf 


AC, 


DW, 


AC, 


sub 


DW 


M„,o ‘' sub, a V 


Then correct C^ at M^ by 


and find 


Ca “ C. - AC^„ 

A A . DW„ 

corr orlg M »o 


-«.T * C„ sina + C. cosa 

DWu Nw ^ A 

M_,o M »a corr 


The total drag is then given by 
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DF 


Sf friction drag determined by the Squire-Young method, 

^ile not perfect, this approach seems to yield good estimates. Also, the 
nvestigation is continuing and the method does not as yet include a correc- 
tion for non-conservative differencing. Finally, ^ used to determine 

the correction factors should be as high as possible.®**® 

Some typical analysis results are shown in figures 7-9, Figure 7 shows 
varying angle of attack at the design Mach number for Airfoil 
109B (critical rooftop design), while figure 8 portrays variations due to 
c anges in M„. As a increases, the flow on the upper surface goes supercrit- 
ical and a shock forms. However, there is a desirable pressure plateau aft 
of the shock wave (ref. 12) which will permit boundary layer recovery. While 
the analysis results indicate no separation for the conditions shown, there 
probably would be shock induced separation at higher a’s. Because of the 
Stratford recovery, such separation would probably lead to a large separated 

region and a very sharp break in the C (a) curve. 

L 

As can be seen in figure 8, as increases a supersonic bubble forms 
and grows and eventually terminates in a shock wave for ^ 0.79. Aft of 
the supersonic zone, the pressures closely follow the original design dis- 
tribution with no separation. These, and other studies, showed that at a C 
of 0,35, drag divergence occurs at M of 0.78. ^ 

ao 

Similar studies were performed for Airfoil 209 (supercritical rooftop 
design), and the angle, of attack variation is shown in figure 9. Note that 
the pressure distribution variation at positive angles of attack is consid- 
erably different from that of Airfoil 109B (fig. 7) in that a shock wave 
fonts immediately and C increases rapidly. Interestingly, the upper surface 
pressure continues to follow the original Stratford recovery aft of the 
shock wave. Also, other studies indicate possible shock induced separation 
for i 0.77 even at zero angle of attack. 

Aaalysis results such as these can also be used to compare airfoils ob- 
tained using different design philosophies. An example for the two designs 
being considered here is shown in figure 10. Notice that for Airfoil 109B, 

Cq is relatively constant up to C,^’g of 0.55, while Airfoil 209 exhibits a’ 
steady increase in C^. in addition, other calculations indicate that 109B 
has little or no drag creep for 0.5 < < 0.7i and C, < 0.5. However, Air- 

209 has 6-18 counts of drag creep in the same range. 

Now it should be pointed out that these results are not critical of 
supercritical airfoils. In fact. Airfoil 109B, whose shape was determined by 
the inverse method using a "critical" rooftop could be used at Cj 's up to 
0.55 with low drags. Thus, it is in essence a supercritical airfoil and 
could be used at design lift coefficients up to 0.55. The point is that with 
the design philosophy and assumptions used in these examples, the inverse 
technique appears to yield the best results by designing the shape with a 
critical rooftop at a lovj The airfoil then can, at least in this case, 
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be used at higher lift coefficients. 


APPLICA'iTON TO THE HIGH LIFT CASE 


technique could be^LccSsfilly’^Ipplled ^rth^r*^^*^ direct-inverse 

specifying the separation point he wa <5 aKi lift case. By 

with experiment by solving thfL^U oLf H 1°. excellent agreement 
boundary conditions upstrfam of separLion and^r 

(pressure specified) downstream of separation Thnr%h°“"‘*®‘'^ conditions 
could similar results be obtained usiL tha ?’n^ ’ *l“®stion arose — 

lfTe\lZToT^ PoLl'f L\e™IneT:s“^ 

in fp: bTe:^p\Xr\1 - ^^own 

surface, the flowfield d 0 «-aY-m 4 j PP opriately modified. On the lower 

specified) including the effects ^of conditions (airfoil 

surface, the flowfield is also computed directly ?" “PP®>^ 

to the separation point, which Is dprat^T^ a viscous Interaction up 

stream of separatiL invlrsfhilf !? solution. Down- 

pressure is Lsumer^o bTconstant J^th'rf f® ®**^ 

have shown that the separated zone pressure which 1 ^"®* studies 

solution, must be determined by conditions at bo P®*-t of the 
at the trailing edge and not iust on rond^f . ? separation point and 

This result is in agreement with the vicinity of separation, 

conditions at the domTtTeL tnd Gross (ref. 14) that 

pressure. For the present^tudi sepaii^tlon bubble determine bubble 

the separated pressure by ^ ^ found adequate to approximate 


- 2 (^ 



sep 

Tit TTTiTizrTT- Ttzx 

with a constant pressure across the Dse^Hn^^^ assumed to be invJscld 

and lower displacement surfaces. formed by the upper 

separation point locatlon'”and^thus”the lift >>ave indicated that the 

boundary layer transition point!^ Thus th! v<=* dependent upon the 

IKAUDES has been modified L ijilul^an 

Puted by a compressible Ihwaites method). naturaJ^tr^L^tLTlL'^Je" i“"' 
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turbulent boundary layer computed by the Nash-Macdonald method. For those 
cases where laminar separation occurs prior to separation, a long or short 
bubble, depending upon local flow conditions, is empirically modeled, and 
then transition is assumed. 

The calculation procedure uses the same iterative successive column re- 
laxation scheme used in the basic program except that the separation point 
and separated pressure level are permitted to vary. A convergence history 
for a typical case is shown in figure 12. Initially some oscillation occurs 
on each grid; but, as can be seen, the values quickly converge. Normally, 
four hundred iterative cycles are performed on both the medium and fine grids. 
The former normally yields 66 points on the airfoil, while the latter yields 
130. 


Results for a GA(W)-2 airfoil are shown in figures 13 and 14. In both 
cases, the lower surface remained entirely laminar, although results with an 
all turbulent lower surface boundary layer showed no significant differences. 
On the upper surface transition with a short separation bubble occurred near 
the leading edge. In general, comparison with experimental data (ref. 18) 
is good with respect to separation point, separated pressure level, and 
C ; and thus the method is quite promising. 

Lm 

Figure 15 shows a comparison with experiment of C versus angle of attack 
for the same airfoil. Similar results have also been obtained at other 
Reynolds numbers. At the present time, research on this approach to the high 
lift problem is continuing in order to ensure that can be determined 

accurately. max 

As a final note, this procedure has also been applied to Airfoil 109B 
discussed previously. Surprisingly, the design Stratford recovery used at 
Mach 0.7A seems to also affect the low speed flow since a sharp break occurred 
in C (a) at 19 degrees and a of 2.08, indicating another disadvantage of 
using a Stratford recovery. 

CONCLUDING REMARKS 


Based upon the results presented here, the following remarks can be 
stated: (1) The present viscous analysis method (TRANDES) is suitable for 

engineering estimates of transonic airfoil data; (2) The present inverse 
design method accounts for the effects of weak viscous interaction in the air- 
foil design process and Is numerically consistent with analysis results; 

(3) The complete potential flow equation coupled with a boundary layer method 
can be used in a direct-inverse fashion to accurately compute the flow about 
airfoils at low speeds having massive separation. 
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i\PPENDIX 


SYMBOLS 


DF 


DW 


M 


'N 


M 

00 

RN 


DW 


AC, 

6* 

Ax 

<P 


axial force coefficient 

two-dimensional drag coefficient 

profile drag coefficient 

wave drag coefficient 

two-dimensional lift coefficient 

two-dimensional quarter chord moment coefficient 

normal force coefficient 

pressure coefficient 

chord 

freestream Mach number 
freestream Reynolds number 
angle of attack 

wave drag coefficient correction 
boundary layer displacement thickness 
trailing edge thickness 
length of separated region 
perturbation potential 


Subscripts 

corr corrected value 

ITE trailing edge 

sep reparation point or region 

sub, a case where flow is entirely subcrltical and at angle of attack a 
orig uncorrected original value 

M ,a at angle of attack rt and Mach number M 

00 ' ” OD 
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Figure 5.- Variation of trailing-edge ordinates with nose-shape thickness 

for upper and lower surfaces. 



Figure 6.- Profile shape and comparison of design and analysis pressure 

distributions for CRAM-209. 
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Figure 7.- Comparison of pressure distributions at four angles 

of attack for CRAM-109-B. 



Figure 8.- Comparison of pressure distributions at 
four Mach numbers for CRAM-109-B. 
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Figure 11 .- Problem formulation. 
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Figure 13.- Theoretical and experimental pressure-distribution 
comparisons. Laminar turbulent case; cx = 12 . 



Figure 14.- Theoretical and experimental pressure-distribution 
comparisons. Laminar turbulent case; a - 18 . 
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Figure 15.- Comparison of theory and experiment for plotted 
against ot for a GA(W)-2 airfoil. 




